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The decomposition of the linear-order metric perturbation is discussed in the context of 
the higher-order gauge-invariant perturbation theory. We show that the linear order metric 
perturbation is decomposed into gauge-invariant and gauge- variant parts on the general back- 
ground spacetime which admits ADM decomposition. This decomposition was an important 
premise of the general framework of the higher order gauge-invariant perturbation theory 
proposed in the papers [K. Nakamura, Prog. Theor. Phys. 110 (2003), 723; ibid. 113 (2005), 
481]. This implies that we can develop the higher-order gauge- invariant perturbation theory 
on generic background spacetime. Remaining issues to complete the general-framework of 
the higher-order gauge-invariant perturbation theories are also discussed. 

§1. Introduction 

Perturbation theories are powerful techniques in many area of physics and the 
developments of perturbation theories lead physically fruitful results and interpreta- 
tions of natural phenomena. 

In physics, physicists want to describe realistic situations in a compact manner. 
Exact solutions in a theory for physical situations are candidates which can describe 
realistic natural phenomena. However, in many theories of physics, realistic situa- 
tions are too complicated and often difficult to describe by an exact solution of a 
theory. This difficulty may be due to the fact that exact solutions only describe 
special cases even if the theory is appropriate to describe the natural phenomena, or 
may be due to the lack of the applicability of the theory itself. Even in the case where 
an exact solution of a theory well describes a physical situation, the properties of the 
physical system will not be completely described only through the exact solution. 
In natural phenomena, there always exist "fluctuations". In this case, perturbative 
treatments of the theory is a powerful tools and physicists investigate perturbative 
approach within the theory to clarify the properties of fluctuations. 

General relativity is a theory in which the construction of exact solutions is not 
so easy. Although there are many exact solutions to the Einstein equation^ 1 these are 
often too idealized. Of course, there are some exact solutions to the Einstein equation 
which well-describe our universe, or gravitational field of stars and black holes. These 
exact solutions by itself do not describe the fluctuations around these exact solutions. 
To describe these fluctuations, we have to consider the perturbations around these 
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exact solutions. Therefore, general relativistic linear perturbation theory is a useful 
technique to investigate the properties of fluctuations around exact solutionsP 

On the other hand, higher-order general-relativistic perturbations also have very 
wide applications. In these applications, second-order cosmological perturbations are 
topical subjeclP^ due to the precise measurements in recent cosmology.® Higher- 
order black hole perturbations are also discussed in some literatureP 1 Moreover, 
as a special example of higher-order perturbation theory, there are researches on 
perturbations of a spherical star^ 1 motivated by the researches on the oscillatory 
behaviors of a rotating neutron star. Thus, there are many physical situations to 
which general relativistic higher-order perturbation theory should be applied. 

As well-known, general relativity is based on the concept of general covariance. 
Intuitively speaking, the principle of general covariance states that there is no pre- 
ferred coordinate system in nature, though the notion of general covariance is math- 
ematically included in the definition of a spacetime manifold in a trivial way. This 
is based on the philosophy that coordinate systems are originally chosen by us, 
and that natural phenomena have nothing to do with our coordinate system. Due 
to this general covariance, the "gauge degree of freedom", which is an unphysical 
degree of freedom of perturbations, arises in general-relativistic perturbations. To 
obtain physically results, we have to fix this gauge degrees of freedom or to extract 
some invariant quantities of perturbations. This situation becomes more compli- 
cated in higher-order perturbation theory. In some linear perturbation theories on 
some background spacetimes, there are so-called gauge-invariant perturbation theo- 
ries. In these theories, one may treat only variables which are independent of gauge 
degree of freedom without any gauge fixing. Therefore, it is worthwhile to investi- 
gate higher-order gauge-invariant perturbation theory from a general point of view 
to avoid gauge issues. 

According to these motivation, the general framework of higher-order general- 
relativistic gauge-invariant perturbation theory has been discussed in some papers^ 1 '^ 1 
by the present author. We refer these works as KN2003 11 ' and KN2005. 12 ' Al- 
though these development of higher-order perturbation theory was originally moti- 
vated by the research on the oscillatory behavior of a self-gravitating Nambu-Goto 
membranepS these works are applicable to cosmological perturbations and we clar- 
ified the gauge-invariance of the second-order perturbations of the Einstein equar 
tions j3JHin3 i n this paper, we refer these works as KN2007® and KN2009JS 

In KN2003p^ we proposed the procedure to find gauge-invariant variables for 
higher-order perturbations on a generic background spacetime. This proposal is 
based on the single assumption that we have already known the procedure to find 
gauge-invariant variables for the linear-order metric perturbation. Under the same 
assumption, we summarize some formulae for the second-order perturbations of the 
curvatures and energy-momentum tensor for the matter fields in KN2009^ and 
KN2009. 14 ' In KN2007P we develop the second-order gauge-invariant cosmological 
perturbation theory after confirming that the above assumption is correct in the 
case of cosmological perturbations. Through these works, we find that our general 
framework of higher-order gauge-invariant perturbation theory is well-defined except 
for the above assumption for linear-order metric perturbations. Therefore, we pro- 
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posed the above assumption as a conjecture in KN2009. 14 ' If this conjecture is true, 
higher-order general-relativistic gauge-invariant perturbation theory is completely 
formulated on generic background spacetime and has very wide applications. 

The main purpose of this paper is to give a proof of this conjecture using the 
premise that the background spacetime admits ADM decomposition. Although some 
special modes are excluded in the proof in this paper, we may say that the above 
conjecture is almost correct for linear-order perturbations on generic background 
spacetimes. This paper is the complete version of our previous short letter 

The organization of this paper is as follows. In £J21 we review the general 
framework of the second-order gauge-invariant perturbation theory developed in 
KN2003™ a nd KN2005 12 -' with some additional explanations. In the context of 
this general framework, the above conjecture is also declared as Conjecture 12.11 in 
this section. In $3l we give a proof of Conjecture 12.11 From pedagogical point of 
view, we consider three different situations of the geometry of the background space- 
time in terms of ADM decomposition. The first situation is trivial ( §3.ip . in which 
we may choose the unit lapse function a = 1, the vanishing shift vector (3 l = 0, and 
the vanishing extrinsic curvature = 0. Through this trivial case, we give the 
essential outline of the proof in more generic situations. The second situation is the 
case where a = 1, f3 l = 0, but ^ ( fl3.2|) . Through this second case, we show 
a technical issue to prove the conjecture 12.11 in terms of ADM decomposition. The 
final situations is most generic case, in which a / 1, ^ ^ 0, Kij 7^ ( §3.3j) . The 
calculations in this case is complicated. However the essential outline of the proof 
is same as in the first trivial case, and the essential technique using in the proof is 
already given in the second case. We also note that we assume that the existence of 
Green functions for two elliptic differential operators in these proofs. The compari- 
son with the proof in the case of cosmological perturbations shown in KN2007" is 
discussed in fgj The final section is devoted to summary and discussion. 

We employ the notation of KN2003 and KN2005 and use abstract index nota- 
tion.^ We also employ natural units in which Newton's gravitational constant is 
denoted by G and the velocity of light satisfies c = 1 . 

§2. General framework of the higher-order gauge-invariant 

perturbation theory 

In this section, we review the general framework of the gauge-invariant perturba- 
tion theory developed in KN2003^D> and KN200dI3 to emphasize that Conjecture O 
is the important premise of our general framework. In §2.1 L we review the basic un- 
derstanding of the gauge degree of freedom in general relativistic perturbation theory 
based on the work of Stewart et alP^ and Bruni et al.® When we consider perturba- 
tions in any theory with general covariance, we have to exclude these gauge degrees 
of freedom in the perturbations. To accomplish this, gauge-invariant variables of 
perturbations are useful, and these are regarded as physically meaningful quantities. 
In §2.2i we review the procedure for finding gauge-invariant variables of perturba- 
tions, which was developed in KN2003. 11 ' After the introduction of gauge-invariant 
variables, in £ j2.3l we review the general issue of the gauge-invariant formulation for 
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the second-order perturbation of the Einstein equation developed in KN2005p^ We 
emphasize that the ingredients of this section do not depend on the details of the 
background spacetime, if the decomposition coni ectur e 12 . 1 1 for the linear-order metric 
perturbation is correct. 

2.1. Gauge degree of freedom in perturbation theory 
2.1.1. Basic idea 

Here, we explain the concept of gauge in general relativistic perturbation theory. 
To explain this, we first point out that, in any perturbation theory, we always treat 
two spacetime manifolds. One is the physical spacetime M, which we attempt 
to describe in terms of perturbations, and the other is the background spacetime 
Mo, which is a fictitious manifold prepared for perturbative analyses by hand. We 
emphasize that these two spacetime manifolds M and Mo are distinct. Let us denote 
the physical spacetime by {M,g a b) and the background spacetime by {Mo, gab), 
where g a b is the metric on M, and g a b is the metric on Mo- Further, we formally 
denote the spacetime metric and the other physical tensor fields on the physical 
spacetime by Q and its background value on the background spacetime by Qo- 

Second, in any perturbation theories, we always write equations for the pertur- 
bation of the physical variable Q in the form 

QCY) = Qo{p) + SQ{p)- (2-1) 

Usually, this equation is simply regarded as a relation between the physical variable 
Q and its background value Qo, or as the definition of the deviation 5Q of the physical 
variable Q from its background value Qo- However, Eq. (|2-ip has deeper implications. 
Keeping in our mind the above fact that we always treat two different spacetimes, 
{M, g a b) and {Mo, 9ab), m perturbation theory, Eq. (|2-1|) is a rather curious equation 
in the following sense: The variable on the left-hand side of Eq. (|2-ip is a variable 
on the physical spacetime {M,g a b), while the variables on the right-hand side of 
Eq. (|2-1|) are variables on the background spacetime, {Mo,g a b)- Hence, Eq. (|2-1|) 
gives a relation between variables on two different manifolds. 

Further, we point out the fact that, through Eq. (|2-ip . we have implicitly iden- 
tified points in two different manifolds {M,g a b) and {Mo, gab)- More specifically, 
Q{ u p") on the left-hand side of Eq. ([24]) is a field on M, and "p" G M. Similarly, 
we should regard the background value Qo{p) of and its deviation SQ{p) of 

Q("p") from Qo{p), which are on the right-hand side of Eq. (|2-ip . as fields on Mo, and 
p € Mo- Because Eq. (|2-ip is regarded as an equation for field variables, it implicitly 
states that the points "p" G M and p E .Mo are same. Therefore, through Eq. (|2-1|) . 
we implicitly assume the existence of a map Mo — > M : p € Mo *— > "p" € M, which 
is called a gauge choice in perturbation theory!^ 

Further, we have to note that the correspondence between points on Mo and 
M, which is established by such a relation as Eq. (|2-ip . is not unique to the per- 
turbation theory with general covariance. Rather, Eq. (|2-1|) involves the degree of 
freedom corresponding to the choice of the map X : Mo \— > M. This is called the 
gauge degree of freedom in general relativistic perturbation theory. Such a degree of 
freedom always exists in perturbations of a theory with general covariance. "Gen- 
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eral covariance" intuitively means that there is no preferred coordinate system in 
the theory. If general covariance is not imposed on the theory, there is a preferred 
coordinate system in our nature, and we naturally introduce this coordinate system 
onto both Mo and M. Then, through this preferred coordinate system, we can 
choose the identification map X. However, due to general covariance, there is no 
such coordinate system in general relativity, and we have no guiding principle to 
choose the identification map X . Actually, we may identify "p" 6 M with q € Mo 
(q p) instead of p G Mq. In the above understanding of the concept of "gauge" in 
general-relativistic perturbation theory, a gauge transformation is simply a change 
of the identification map X. 

These are the basic ideas necessary to understand gauge degree of freedom in the 
general relativistic perturbation theory proposed by Stewart and WalkerPD 1 This 
understanding has been developed by Bruni et al.J^ and by the present author.^'^ 

2.1.2. Formulation of perturbations 

To formulate the above understanding in more detail, we introduce an infinites- 
imal parameter A for the perturbation. Further, we consider the (n + 1) + 1- 
dimensional manifold Af = MxM, where n+1 = dim.M and AGl. The background 
spacetime Mo = A~\ x=0 and the physical spacetime M = M\ = Af\ R=x are also 
submanifolds embedded in the extended manifold Af. Each point on Af is identified 
by a pair, (p, A), where p G M\, and each point in the background spacetime A4o in 
Af is identified by A = 0. 

Through this construction, the manifold Af is foliated by (n + l)-dimensional 
submanifolds M\ of each A, and these are diffeomorphic to the physical spacetime 
M and the background spacetime Mq. The manifold has a natural differen- 
tiable structure consisting of the direct product of M and M. Further, the perturbed 
spacetimes M\ for each A must have the same differential structure with this con- 
struction. In other words, we require that perturbations be continuous in the sense 
that (M,g a b) and (Mo, gab) are connected by a continuous curve within the ex- 
tended manifold Af. Hence, the changes of the differential structure resulting from 
the perturbation, for example the formation of singularities, are excluded from our 
consideration. 

Let us consider the set of field equations 

£[Qx] = (2-2) 

on the physical spacetime M\ for the physical variables Q\ on M\. The field 
equation (|2-2p formally represents the Einstein equation for the metric on M\ and 
the equations for matter fields on M\. If a tensor field Q\ is given on each M\, 
Qx is automatically extended to a tensor field on N by Q(p, A) := Q\(p), where 
p € M\. In this extension, the field equation (|2-2p is regarded as an equation on 
the extended manifold Af. Thus, we have extended an arbitrary tensor field and the 
field equations (|2-2j) on each M\ to those on the extended manifold M. 

Tensor fields on Af obtained through the above construction are necessarily "tan- 
gent" to each M\, i.e., their normal component to each M\ identically vanishes. To 
consider the basis of the tangent space of Af, we introduce the normal form and 
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its dual, which are normal to each A4\ in J\f. These are denoted by (d\) a and 
(d/d\) a , respectively, and they satisfy (dX) a (d/dX) a = 1. The form (dX) a and its 
dual, (d/d\) a , are normal to any tensor field extended from the tangent space on 
each A4\ through the above construction. The set consisting of (dX) a , (d/dX) a , and 
the basis of the tangent space on each A4\ is regarded as the basis of the tangent 
space of Af. 

To define the perturbation of an arbitrary tensor field Q, we compare Q on the 
physical spacetime A4\ with Qq on the background spacetime, and it is necessary 
to identify the points of Ad\ with those of A4q. This point identification map is the 
so-called gauge choice in the context of perturbation theories, as mentioned above. 
The gauge choice is made by assigning a diffeomorphism X\ : Af — > TV such that X\ 
: A4q — > A4\. Following the paper of Bruni et al.p we introduce a gauge choice 
X\ as an exponential map on Af, for simplicity. We denote the generator of this 



exponential map by x rf. 



This generator x rj a is decomposed by the basis on the 



tangent space of Af which are constructed above. The arbitrariness of the gauge 
choice X\ is represented by the tangential component of x rf to AA\. 

The pull-back X^Q, which is induced by the exponential map X\, maps a tensor 



field Q on to a tensor field X^Q on Mo. In terms of this generator r] a , the 
pull-back X? Q is represented by the Taylor expansion 



Q(r) = Q{X x {p)) = X* x Q(p) 



Q(p) + \£ Xri Q\ p + -X 



2 £ 2 



+ 0(A 3 ), (2-3) 



where r = X\(p) G A4\. Because p G A4q, we may regard the equation 
X*xQ(p) = Qo(p) + A £ Xv Q\ u „ ( P ) + 7: A 2 £l n Q tA ( P ) + C»(A 



Mo 



(2-4) 



as an equation on the background spacetime A4q, where Qq = Q\j^ is the back- 
ground value of the physical variable of Q. Once the definition of the pull-back of 
the gauge choice X\ is given, the perturbations of a tensor field Q under the gauge 
choice X\ are simply defined by the evaluation of the expansion (|2-4p on A^o 



X{Q 



^\Mo 



Qo + A ( i } Q + ^A 2 ^Q + 0(A 3 



1 



2(2), 



(2-5) 



i.e., 



Mo 



y<5 := 



£ 2 



XjJ Q 



Mo 



(2-6) 



We note that all variables in this definition are defined on .Mo- 

2.1.3. Gauge transformation 

Here, we consider two different gauge choices. Suppose that X\ and 3^a are two 
exponential maps with the generators x r] a and ^if on M, respectively. The integral 
curves of each x rf and ^rf in M are the orbits of the actions of the gauge choices 
X\ and 3^A) respectively. Since we choose x r\ a and ^yf so that these are transverse to 
each A4\ everywhere on Af, the integral curves of these vector fields intersect with 
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each A4\. Therefore, points lying on the same integral curve of either of the two are 
to be regarded as the same point within the respective gauges. When these curves 
are not identical, i.e., the tangential components of x rj a and ^rf to each A4\ are 
different, these point identification maps X\ and y x are regarded as two different 
gauge choices. When we have two different gauge choice X\ and y x , we have two 
different representations of the perturbative expansion of the pulled-backed variables 
X*Q x \ Mo and y* x Qx\ Mo : 



y*xQx\ Mn = Qo + *yQ + l^yQ + 0(A 3 ), (2-8) 



Although these two representations of the perturbations are different from each other, 
these should be equivalent because of general covariance. 

Now, we consider the gauge-transformation rules between two different gauge 
choices. In general, the representation x Q\ on of the perturbed variable Q on 
A4\ depends on the gauge choice X x . If we employ a different gauge choice, the 
representation of Q\ on M.q may change. Suppose that X x and y x are two different 
gauge choices and the generators of these gauge choices are given by x r\ a and -^r] a , 
respectively. In this situation, the change of the gauge choice from X\ to y x is 
represented by the diffeomorphism 

$\:={XxY X °y\- (2-9) 

This diffeomorphism <P X is the map ^\ : A4q — > A4q for each value of A £ M. The 
diffeomorphism <P\ does change the point identification, as expected from the under- 
standing of the gauge choice discussed above. Therefore, the diffeomorphism <P X is 
regarded as the gauge transformation <P X : X\ — > y x . 

The gauge transformation <P\ induces a pull-back from the representation X Q\ 
of the perturbed tensor field Q in the gauge choice X\ to the representation ^Q x in 
the gauge choice y\. Actually, the tensor fields X Q\ and ^Q\, which are defined on 
M.q, are connected by the linear map <P* X as 



y Qx = yiQ\ Ma = (yiixxx-yq 



Mo 



= {X^yxY (X* X Q) = <P* X Qx- (2-10) 

Mo 

According to generic arguments concerning the Taylor expansion of the pull-back of 
a tensor field on the same manifoldf^® it should be expressed the gauge transfor- 
mation <P x x Qx m the form 

@*x X Q = X Q + A%*Q + y {% + 4 } X Q + 0(A 3 ), (2-11) 

where the vector fields and are the generators of the gauge transformation <P X . 

Comparing the representation (|2- and that in terms of the generators x rj a and 
^r] a of the pull-back y^ o (X^ ls ) X Q (= <P* X X Q), we obtain explicit correspondence 
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y 

(2), 



between and { x n a , y i] a } as follows: 

£? = V-V, e 2 a = [^,^] a . (2-i2) 

Further, because the gauge transformation <P\ is a map within the background space- 
time Mo, the generator should consist of vector fields on Mo. 

We can now derive the relation between the perturbations in the two different 
gauges. Up to second order, these relations are derived by substituting (|2-7p and 
(1771) into PTT]> : 

Mq_Mq = £ 5i q , (2-13) 

'Q - $Q = 2£ ?(1) «Q + {£ k2) + Qo- (2-14) 

Here, we comment on the generic formula for the Taylor expansion (|2-11|) . In the 
case where we regard the pull-backs X£ and y% of the gauge choices are exponential 
maps, the product of two exponential maps is also written by the exponential of the 
infinite sum of the Lie derivative along infinitely many generators through Baker- 
Campbell-Hausdorff formulaP^'^ These infinitely many generators are constructed 
by the commutators of x rf and ^r] a , which are regarded as higher-order derivatives 
of x rf and ^rf in TV and are regarded as the vector fields on M . The expression 
(|2Tip of the Taylor expansion is just the expression up to 0(A 3 ) of the Baker- 
Campbell-Hausdorff formula. Of course, we may generalize the gauge choice XX 
and y\ to more general class of diffeomorphism than the exponential map. Even 
in this case, the Taylor expansion (|2Tip is correct. Since the gauge-transformation 
rules (|2-13p and (|2T4p are direct consequences of the Taylor expansion (|2Tip , these 
gauge-transformation rules are not changed even if we generalize the gauge choice 
X£ and and we may regard that two generator and £? 2 ^ in Eqs. (|2T3p are 
independent of each other. Therefore, we may say that gauge-transformation rules 
(|2-13p and (|2-14|) are most general gauge-transformation rules of the first and second 
order, respectively. 

2.1.4. Gauge invariance 

We next introduce the concept of gauge invariance. The gauge invariance con- 
sidered in this paper is order by order gauge invariance proposed in KN2009P* We 

(p) 

call the /cth-order perturbation X Q is gauge invariant iff 

{ 5Q = {k y ] Q (2-15) 

for any gauge choice X\ and y\. Through this concept of order by order gauge 
invariance, we can decompose any perturbation of Q into the gauge-invariant and 
gauge-variant parts, as shown in KN2003.^ In terms of these gauge-invariant vari- 
ables, we can develop the gauge-invariant perturbation theory. However, this devel- 
opment is based on a non-trivial conjecture, i.e., Conjecture 12.11 for the linear order 
metric perturbation as explained below. 

2.2. Gauge-invariant variables 

Inspecting the gauge-transformation rules (|2-13j) and (|2T4p . we define gauge- 
invariant variables for metric perturbations and for arbitrary matter fields. First, we 
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consider the metric perturbation and expand the metric g a b on M, which is pulled 
back to M.q using a gauge choice X\ in the form given in (|2-5p . 

A 2 

X\9ab = 9ab + Xrflab + y dab + 0'\\), (2-16) 

where g a b is the metric on A4q. Of course, the expansion (|2-16|) of the metric depends 
entirely on the gauge choice X\. Nevertheless, henceforth, we do not explicitly 
express the index of the gauge choice X\ if there is no possibility of confusion. 

Our starting point to construct gauge-invariant variables is the following conjec- 
ture^ for the linear-order metric perturbation h a b defined by Eq. (|2-16p : 
Conjecture 2.1. If there is a tensor field h a b of the second rank, whose gauge 
transformation rule is 

yhab — xhab = ££ {1) 9ab, (2"17) 

then there exist a tensor field Hab and a vector field X a such that h a b is decomposed 
as 

h a b =■ Hab + £x9ab, (2" 18) 

where % a b and X a are transformed as 

yriab - xHab = 0, yX a - x X a = £^ (2-19) 

under the gauge transformation \2-13\) , respectively. 

In this conjecture, H a b is gauge-invariant in the sense as mentioned above, and 
we call H a b as gauge-invariant part of the linear-order metric perturbation h a b- On 
the other hand, the vector field X a in Eq. (|2-18|) is gauge dependent, and we call X a 
as gauge-variant part of the metric perturbation h a b. 

The main purpose of this paper is to prove Con j ectur e 1 2 . 1 1 in some sense. In the 
case of the cosmological perturbations on a homogeneous and isotropic universe, we 
confirmed Conjecture 12.11 is correct except for some special modes of perturbations, 
and then we developed the second-order cosmological perturbation theory in a gauge- 
invariant manner.^ On the other hand, in the case of the perturbation theory on 
a generic background spacetime, this conjecture was highly non-trivial due to the 
non-trivial curvature of the background spacetime. We see this situation in detail in 
^31 However, before going to the proof of Conjecture 12. 1[ we explain how the higher- 
order gauge-invariant perturbation theory is developed based on this conjecture, 
here. Through this explanation, we emphasize the importance of Conjecture 12. 1[ 

As shown in KN2003p* the second-order metric perturbations l a b are decom- 
posed as 

Lb ='■ £ab + 2£xh a b + [£y - £\) 9ab, (2-20) 

where C a b and Y a are the gauge-invariant and gauge-variant parts of the second 
order metric perturbations, i.e., 



y£ab - x£ab = 0, yY a - X Y a = £" 2 ) + X] 



(2-21) 
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Actually, using the gauge- variant part X a of the linear-order metric perturbation 
h a b, we consider the tensor field L a f, defined by 

L a b ■= lab — 2£ x h ao + £ 2 x9ab- (2-22) 

Through the gauge-transformation rules (|2-14p and (|2-19p for l a b and X a , respec- 
tively, the gauge-transformation rule for this variable L a b is given by 

yL a b - xL a b = £a9ab, 0"° '■= £(2) + [£(1)>^] ° • (2-23) 

This is identical to the gauge-transformation rule (|2-17p in Conjecture 12.11 and we 
may apply Conjecture 12.11 to the variable L ao . Then, L a b can be decomposed as 

L a b = £ab + £Y9ab, (2-24) 

where the gauge-transformation rules for C a b an d Y a are given by Eqs. (|2-2ip . To- 
gether with the definition (|2-22p of the variable L a \,i the decomposition (|2-24|) leads 
the decomposition (|2-20p for the second-order metric perturbation l a b- 

Furthermore, as shown in KN2003, using the first- and second-order gauge- 
variant parts, X a and Y a , of the metric perturbations, the gauge-invariant variables 
for an arbitrary tensor field Q other than the metric are given by 

(1) Q := - £xQo, (2-25) 

(2) Q : = ( 2 )Q - 2£ X ^Q -{£ y - £ 2 x) Qq. (2-26) 

It is straightforward to confirm that the variables *P'Q defined by ()2-25p and ()2-26p are 
gauge invariant under the gauge-transformation rules f)2- 13j) and ()2-14j) . respectively. 

Equations (|2-25|) and (|2-26|) have an important implication. To see this, we 
represent these equations as 

«Q = WQ + £ X Q , (2-27) 

(2) Q = (2) Q + 2£ X ^Q + {£ Y - £ 2 x}Qo- (2-28) 

These equations imply that any perturbation of first and second order can always be 
decomposed into gauge-invariant and gauge- variant parts as Eqs. (|2-27p and (!2-28p . 
respectively. 

2.3. Second-order gauge-invariant perturbation theory 

When we consider the first- and the second-order perturbations of the Einstein 
equation, we have to consider the perturbative expansion of the Einstein tensor and 
the energy momentum tensor. Now, we consider the perturbative expansion of the 
Einstein tensor on A4 \ as 

G a b = G a b + A«G a 6 + \\^ 2) G a b + 0(A 3 ). (2-29) 

As shown in KN2005, 12 ' the first- and the second-order perturbation of the Einstein 
tensor are given by 

(1) G a h = ^G a b [H\ + £xG a \ (2-30) 

(2) G a b = [C] + ^g a b [H,H] + 2£ x ^G b + {£ Y - £\) G a b , (2-31) 
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where 



{1) Ga b [A] 

{2) G a b [AB] 
M£ a b [A,B] 



Vz a b [A]- l ~5 a h Vz c c [A], ( 

-2S/ [a H d] bd [A]-A cb R ac , ( 

^E a b [A,B]-ls a b ^E c c [A,B], ( 

2R ad B c ( b A d ^ + 2H [a d * [A] H d] b e [B] + 2H [a de [B] H d b e [A] 
+2A e d V [a H d] be [B] + 2B e d V [a H d] be [A] 
+2A b V [a H d] cd [B] + 2B c b V [a H d] cd [A] , ( 



(2-32) 



(2-33) 



(2-34) 



(2-35) 



and 



H ab c [A] 

H abc [A] 
H a \ [A] 



9cd H ab d [A], H a bc [A]:=g bd H ad c [A] 
9cdH a hd [A]. 



(2-36) 



(2-37) 



We note that W G a b [*] and {2) G a b {*,*\ in Eqs. (12^301) and (12^311) are the gauge- 
invariant parts of the perturbative Einstein tensors, and Eqs. (|2-30p and ()2-3ip have 
the same forms as Eqs. (|2-27j) and (|2-28p . respectively. 

We also note that ^Q a b [*] and ^Q b [*, *] defined by Eqs. (f2^32V(l2^35D satisfy 
the identities 



for arbitrary tensor fields A ab and B ab , respectively. We can directly confirm these 
identities without specifying arbitrary tensors A ab and B ab of the second rank, re- 
spectively. These identities (|2 38|) and (|2-39p guarantee the first- and second-order 
perturbations of the Bianchi identity V b G b = 0. This implies that our general 
framework of the second-order gauge- invariant perturbation theory is self-consistent. 
On the other hand, the energy momentum tensor on A^^ is also expanded as 



According to Eqs. (|2-27|) and ([2-280 . we can also decompose the first- and the second- 



Va {1) G b a [A] 

V a Wg b a [A,B] 




(2-38) 



(2-39) 



T a b = T a b + X^T a 6 + \\^ 2) T a b + 0(A 3 ). 



(2-40) 




as 



(2-41) 
(2-42) 
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These decompositions are confirmed in the case of a perfect fluid, an imperfect fluid, 
and a scalar field in KN2009PJ Furthermore, in KN2009JS we also showed that 
equations of motion for the matter field, which are derived from the divergence of 
the energy-momentum tensors, are also decomposed into gauge-invariant and gauge- 
variant parts as Eqs. (|2-27p and f|2-28j) . Therefore, we may say that the decomposition 
formulae (|2-27p and (I2-28P are universal. 

Imposing order by order Einstein equations 

G a h = MT a \ ^G a b = 8^>T a b , ^G a b = 8^ 2 >T a b , (2-43) 

the first- and the second-order perturbation of the Einstein equations are automati- 
cally given in gauge-invariant form as 

Wg a b [H] = 8irG^T a \ {1) g a b [C] + ^g a b [H, H] = 8ttG ^T a b . (2-44) 

Furthermore, in KN2009,^ we also showed that the equations of motion for matter 
fields, are automatically given in gauge-invariant form. Thus, we may say that any 
equation of order by order is automatically gauge-invariant and we do not have to 
consider the gauge degree of freedom at least in the level where we concentrate only 
on the equations of the general relativistic system. 

We can also expect that the similar structure of equations of the systems will be 
maintained in the any order perturbations and our general framework be applicable 
to any order general-relativistic perturbations. Actually, decomposition formulae for 
the third-order perturbations in two-parameter case which correspond to Eqs. (|2-27j) 
and (|2-28p are given in KN2003.^ Therefore, similar development is possible for the 
third-order perturbations. Since we could not find any difficulties to extend higher- 
order perturbations^ except for the necessity of long cumbersome calculations, we 
can construct any order perturbation theory in gauge-invariant manner, recursively. 

We have to emphasize that the above general framework of the higher-order 
gauge-invariant perturbation theory are independent of the explicit form of the back- 
ground metric g a b, except for Conjecture 12. 1[ and are valid not only in cosmological 
perturbation case but also the other generic situations if Conjecture 12. II is true. This 
implies that if we prove Conjecture 12.11 for the generic background spacetime, the 
above general framework is applicable to perturbation theories on any background 
spacetime. This is the reason why we proposed Conjecture 12.11 in KN2009.^ 

Thus, Conjecture 12.11 is the important premise of our general framework of 
higher-order gauge-invariant perturbation theory. In the next section, we give a 
proof of Conjecture 12.11 on the generic background spacetime which admits ADM 
decomposition (see Appendix |A|) . 

§3. Decomposition of the linear-order metric perturbation 

Now, we give a proof of Conjecture 12. II on general background spacetimes which 
admit ADM decomposition (see Appendix [A]). Therefore, the background spacetime 
A4o considered here is n + 1-dimensional spacetime which is described by the direct 
product M x U. Here, R is a time direction and U is the spacelike hypersurface with 
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dim 17 = n embedded in Mq. This means that Mq is foliated by the one-parameter 
family of spacelike hypersurface £{t), t G M is a time function. The metric on A4q 
is given as Eq. (|A-23p . i.e., 

9ab = -a 2 {alt) a (dt) h + q ij (dx i + p i dt) a (dx j + ftdt) b , (3-1) 

where a is the lapse function, /3* is the shift vector, and q ab = qij(dx l ) a (dx l ) b is the 
metric on £(t). The inverse of Eq. (|3-1|) is given by Eq. (|A-2ip in the Appendix lAl 
To consider the decomposition (|2-18p of h ab , first, we consider the components 
of the metric h ab as 

Kb = h tt {dt) a (dt) b + 2h ti {dt) {a (dx i ) b) + hijidx^aidx^b. (3-2) 

The components hu, hu, and hij are regarded as a scalar function, a vector field, 
and a tensor field on the spacelike hypersurface S(t), respectively. From the gauge- 
transformation rule (|2-17p . the components {h tt , hu, h^} are transformed as 

yha - xhtt = 2dt£t ~ - {dta + pD ia - &pK io ) & 

+a 2 D l a - ap h D l p k - ppDjo) &, (3-3) 
2 

yh u - xhi = dtii + A& (Diet - P> d K i3 \ £ t 

a 

-| {-a 2 K\ + P j p k K ki - fiDia + aDrfi) (3-4) 
2 2 k 

yhij — xhij = 2-D(j£j) + —Kij^ t — —j3 Ky^, (3-5) 

where is the extrinsic curvature of £ defined Eq. (|A-36p in Appendix [X] and Di 
is the covariant derivative associate with the metric qij (D^qj^ = 0). 

Apparently, the gauge-transformation rules (|3-3p - (|3-5p have the complicated 
form and it seems difficult to find the decomposition as (|2-18p . Therefore, we consider 
the proof of the decomposition (|2-18|) from the simpler situations of the background 
spacetime. We consider the three cases. In H3.1\ the case where a = 1, f3 l , and 
Kij = is considered. In this simplest case, we may consider the non-trivial in- 
trinsic curvature of E. From the view point of a proof of Conjecture 12.11 this case 
is trivial. However, we find the outline of the proof of Conjecture 12.11 through this 
case. Second, in §3.21 we consider the case where a = 1, j3 l = 0, but Kij ^ 0. 
This case includes not only many homogeneous background spacetimes but also the 
Schwarzschild spacetime. Furthermore, we note that the most non-trivial technical 
part of the proof of Conjecture 12.11 is given in this case. Finally, in N3.31 we consider 
the most general case where a / 1, /3 1 / 0, and Km ^ for completion. 

3.1. a = 1, 0i = 0, and Ky = case 

Here, we consider .Mo satisfies the conditions 

a = l, Pi = Q, K ij= 0. (3-6) 
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In this case, the gauge-transformation rules (|3-3p ~ (|3-5p are given by 

yhtt - xht = 2dt£ t , (3-7) 
yhu - xh u = dtii + Aft, (3-8) 
yhij - xhij = 2D(i£jy (3-9) 

To prove Conjecture l2.ll we consider the decomposition of the symmetric tensor field 
on U reviewed in Appendix [Bj 

hti = Afyvx) + h(y)i, D l h {v)i = 0, (3-10) 

hij = -Qij h (L) + h (T)ij, Q Vh (T)ij = 0, (3-11) 

h(T)ij = (Lh(rv))ij + h (TT)ij, D % h(TT)ij = 0, (3-12) 

h(TV)i = D ih(TVL) + h(rw)i, D % h( T VV)i = 0, (3-13) 

where (Lh(xv))ij is defined by [see Eq. (|B- 5|) in Appendix IB] 

2 

(Lh( T v))ij '■= D ih(rv)j + Djh(rpv)i - —qijD h( T v)l- (3-14) 

To derive gauge-transformation rules for {h u , fyyn, fyy)i, fyi), ^(TV)i, ^(TT)y}, or 
{^tt, /i(VL), ^(y)i, fyz,), hi^vL), h(rvv)h h( T T)ij}, we decompose the generator & as 

£i = Af(£) + f(V)i, £>V)< = 0. (3-15) 

Through the decomposition (I3-I0p and f|3- 15|) . the gauge-transformation rule 
(|3-8|) is given by 

yhu - xhu = Di (yh(y L ) - xhfyV)) + yh(y^ - xh(y)i 

= A ft) +^00« ( 3 - 16 ) 

where we used dtDif = Didtf for an arbitrary scalar function /. Taking the diver- 
gence of Eq. (|3T6p , we see that 

4 (yVl) " AT Vl)) = 4 (3tf(L) + ft) + (3-17) 

Since Kij = 0, a = 1, and = in our case, Eq. (|A-37j> yields 

K l3 = -^d tgij = 0. (3-18) 
From Eqs. (bA35l) . (pMSD . and the property A£ (v)i = 0, we can easily verify that 

D%Z(y)i = d t D% v)i = 0. (3-19) 
Then the gauge-transformation rule ()3-17j) is given by 

A (yh { VL) ~ xh (VL) ) = A (dtt (L) + f t ) . (3-20) 
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Here, we assume the existence of the Green function of the Laplacian A := D l Di 
and ignore the mode which belongs the kernel of the Laplacian A. Then, we obtain 

yh( V L) ~ xh(yL) = dti(L) + it- (3-21) 
Substituting Eq. (|3-2ip into Eq. (|3T6p . we obtain 

yh(y)i - xh(y)i = d t S,(v)i- (3-22) 

Next, through Eq. (|3-11|) . the trace part of Eq. (|3-9|) is given by 

q ij yhij - q^xhij = yh {L) - x h {L) = 2D% = 2A£ (L) , (3-23) 

where we used Eq. (|3-15p . Then, through the decomposition (j3- 12j) . the gauge trans- 
formation rule of the traceless part ^(T)ij is given by 

y h (T)ij ~ xh(T)ij = y(Lh( TV }).. - x (Lh^ TV ^).. + yh( T T)ij ~ xh(TT)ij 

= • (3-24) 
Taking the divergence of Eq. (|3-24p . we obtain 

V ]l (yh (TV)l - x h (TV)l - &) = 0, (3-25) 
where the derivative operator T> lJ is defined by 

V ij := q ij A + M - ~j D i D j + R ij , (3-26) 

and its properties are discussed in Appendix[Bj Here, we assume the existence of the 
Green function of the derivative operator T> % ^ and ignore the modes which belong to 
the kernel of the derivative operator . Then, we obtain 

yh(rv)l - xh(rv)i = (3-27) 
Substituting Eq. <^27\i into Eq. S^M\ . we obtain 

yh(TT)ij ~ xh( TT )ij = 0. (3-28) 
Moreover, through the decomposition (|3- 13|) . we obtain 

yh(TVL) - xh(TVL) = (3-29) 
yh(TVV)i - xh(TVV)i = C(V)i- (3-30) 

Here, we have also ignored the mode which belongs the kernel of A. 

In summary, we have obtained the gauge-transformation rule of variables {hu, 

h(TT)ij} as 

yhtt - xht = 2d t £t, (3-31) 
yh(VL) - xh(y L ) = 9 t C(L) + 6; (3-32) 
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yh(v)i - x\v)i = <h€(y)i> (3-33) 
yh {L) - xh {L) = 2D% = 2A£ {L) , (3-34) 

y h {TVL) ~ xh(TVL) = £(£)> (3-35) 
yh(TVV)i - xh(TVV)i = €(y)h (3-36) 
yh(TT)ij - xh(TT)ij = 0- (3-37) 

Now, we construct gauge-invariant variables. First, the gauge-transformation 
rule (|3-37p shows hi^T)^ ls g au g e invariant by itself: 

Xij '■= h(TT)ij- (3-38) 

From Eq. (j3- 12|) . this variable Xij satisfy the transverse-traceless condition D l Xij = 
= q ij Xi r 

Second, from the gauge-transformation rules (|3-33p and (|3-36j) . we can define the 
gauge-invariant variable for the vector mode as 

Vi := h(y)i - d t h(rvv)i- (3-39) 
Actually, gauge-transformation rule for the variable z/j is given by 

y^i - xn = {yh(v)i - dtyh(rvv)i) ~ [xh(y)i ~ dtxh( T vv)i) 

= dt^ v) i - d t ^ v)i = 0. (3-40) 

From the definition, the gauge-invariant variable i/j satisfy the transverse condition 
D % Vi = as a result. 

Next, we consider scalar modes. First, from the gauge-transformation rules 
(|3-34j) and ()3-35|) . we see that the variable W defined by 

- 2nV := h (L) - 2Ah (TVL) (3-41) 

is gauge invariant. Actually, the gauge-transformation rule for & is given by 

y &- x &= -— ( y h [L) - 2A y h {TVL) ) + — ( x h (L) - 2A x h (TVL) ) = 0.(3-42) 

To define another gauge-invariant variables for scalar modes, we consider the gauge- 
transformation rule (|3-32p and (|3-35p . We find that the variable Xt defined by 

Xt ■= h(yL) ~ d t h( TVL ) (3-43) 

is transformed as 

y X t - x X t = ^. (3-44) 

Using this variable Xt, we define the gauge-invariant combination as 

- 2$ := h tt - 2d t X t . (3-45) 

Actually, we can easily check the variable ^ is gauge invariant under the gauge- 
transformation rules (|3-3ip and (|3-44[) . 
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In terms of the gauge-invariant variables defined by Eqs. (j3-38|) . ()3-39p . (|3-4ip . 
and (|3-45p . and the variable X t defined by Eq. (|3-43p . the original set {htt, hu, hij} 
of the components of the linear metric perturbation is given by 

htt = -2$ + 2d t X t , (3-46) 
hi = n + d t h(rv)i + DiX t , , (3-47) 
hij = -2qij& + Xij + Dih(TV)j + D j h (TV)i- (3-48) 

On the other hand, the components of Eq. (|2T8p with 

X a =:X t (dt) a +Xi{dx% (3-49) 

are given by 

h tt = n tt + 2d t X t , (3-50) 

hti=nti + d t X i + D i Xt, (3-51) 

hij = Uij + DiXj + DjXi. (3-52) 

Since the variable Hu , Hu , and Hi j are gauge- invariant and the gauge-transformation 
rules for the variable Xt and Xi are given by 

yX t - xX t = it, yXi - xXi = (3-53) 

respectively, we may naturally identify the variables as follows: 

Hft = -2#, Hti = Vi, Uij 2<Hj<1' ■ Xij,, ( 3 - 54 ) 
X t = X t , X l = h (TV)l . (3-55) 

The gauge-transformation rules (|3-27p [or equivalently Eqs. f)3-35f) and ()3-36p ] and 
(|3-44p support these identifications. Thus, Eqs. (|3 46|) — f[3-48j) show that the linear- 
order metric perturbation h a b is decomposed into gauge-invariant and gauge- variant 
parts as Eq. (|2T8p in the case for the background spacetime with a = 1, /3* = 0, and 
Kij = 0. 

3.2. a = 1, ft = 0, but ^ case 

Here, we also consider the case where the background metric is described by 
Eq. (|3-ip with a = 1, ft = 0, but the background spacetime has the non-trivial 
extrinsic curvature Kij ^ 0. In this case, the extrinsic curvature K^ is proportional 
to the time derivative of the metric qij on £ 

K ij = - \ d tQij (3-56) 

from Eq. (|A-37p in Appendix [A] and gauge-transformation rules ()3-3p - (j3-5p for the 
components {hu, hu, h^} are given by 

yh tt - xhtt = 2d t it, (3-57) 
yh u - xh ti = d t Zi + D& + 2/v%, (3-58) 
yhj - x hij = 2D (i Z f) + 2K^ t - (3-59) 
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Inspecting gauge-transformation rules (|3-58p - (|3-59p . we first define a new sym- 
metric tensor field H a b whose components are defined by 

Htt '■= hu, H t i := h t i, Hij := — 2KijX t . (3-60) 

Here, we assume the existence of the variable Xt whose gauge-transformation rule is 
given by 

y X t - x X t = (3-61) 

The existence of the variable Xt is confirmed later soon. Similar technique is given 
by T. S. Pereira et al.^ in the perturbations on Bianchi type I cosmology. Now, the 
components Hu and are regarded as a vector and a symmetric tensor on U(t), 
respectively. Then, we may apply the decomposition in Appendix [B] to Hu and H; L y. 

H ti = Dih {VL) + h (v)i , £P/» (v)i = 0, (3-62) 

Hij = ^Qij h (L) + h(T)ij, Q tj h(X)ij = °> ( 3 ' 63 ) 
h (T)ij = [Lh^xv))^ + h {TT)ij, D l h( TT ^j = 0, (3-64) 
h(TV)i = Dih(TVL) + h(rw)h D l h(rvv)i = °- (3-65) 
Through Eqs. (|3 62|) and (|3-6U|) . the gauge-transformation rule (|3-58j) is given by 

A (yh(yL) ~ xh(yL)) + yh(y)i ~ x\v)i = + D& + 2K\^. (3-66) 
Taking the divergence of this gauge-transformation rule (|3-66p . we obtain 

4 (yh iV L) ~ xh(y L) ) = D'dtti + A£ t + 2D* , (3-67) 

where we used the divergenceless property of the variable hry^. 

Now, we evaluate the term D l dt£i in Eq. (|3-67p . First, we note that the time- 
derivative of the intrinsic metric dtq%j is given by Eq. (|3-56p . i.e., 

d tqij = -2Kij. (3-68) 

Keeping the relation (|3-68p in our mind, we can evaluate Didt£j as 

Didtij = dtDitj - C k (DiKj k + DjKi k - D k K^ . (3-69) 

From Eq. (|3-69|) . we further evaluate the term D l dt£,i as 

D'dti, = d t [D%) + ( k D k K - 2Di ■ (3-70) 

Through Eq. (|3-70|1 . the gauge-transformation rule (|3-67p is given by 

^ {yh {V L) ~ xh {VL) ) = d t (D%) + + Z k D k K. (3-71) 
Here again, it is convenient to introduce the decomposition of £j as Eq. (|3T5p . i.e., 

6 = A£ ( l) + £(v> D% v)i = 0. (3-72) 
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Through Eq. ()3-70p . the gauge-transformation rule ()3-7ip is given by 

^ {yh(VL) ~ x\vl) - it) = Ad t i {L) + 2 A (K l W^ {L) ) + D k K^ [v)k . (3-73) 

Since we ignore the modes which belong to the kernel of the operator A, we obtain 
the gauge-transformation rule for the variable fyyn as 

yhvL) ~ xh {VL) = & + dti(L) + A- 1 [2A {K ij D^ {L) ) + D k K£ {v)k ] . (3-74) 

Substituting Eq. (|3-74|) into Eq. (|3-66|) . we obtain the gauge-transformation rule for 
the variable hry\i as follows: 

yh( V )i - xh(y)i = d t C(v)i + 2KJ i D j((L) + 2K] S{v)j 

-DiA- 1 [2A {K ij D 3 i {L) ) + D k K£ iv)k \ . (3-75) 

The divergenceless property of Eq. (|3-75p can be easily checked through Eq. (|3-70p . 
From Eqs. (|3-59p and (|3-6ip . the gauge transformation rule for Hij is given as 

yHij — x&ij = (yhij — 2KijyX t ^j - (xh>ij — 2ify^tJ = ZD^jy (3-76) 

In terms of the decomposition (|3-63p - (|3-65p . the gauge-transformation rules for the 
variables fyn and hrj>)ij is derived from 

—Qij {yh(L) - xh(L)) + {yh(T)ij - xh(r)ij) = 2D {iij)- (3-77) 
The trace part of the gauge transformation rule (|3-77p is given by 

yh iL ) ~ xh {L) = 2D%, (3-78) 
and the traceless part of the gauge transformation rule (|3-77p is given by 

y h (T)ij ~ xh(T)ij = ( L Qij • (3-79) 
Applying Eq. ()3-64p . the gauge-transformation rules ()3-79p are given by 

(L (yh (TV) - xh^ TV ))).. + (yhfTj^ij - xh(TT)ij) = ( L Qij ■ (3-80) 
Taking the divergence of Eq. (|3-80p . we obtain 

V jl [yh (TV)l - x h (TV)l - 6] = 0. (3-81) 

Since we ignore the modes which belong to the kernel of the elliptic derivative oper- 
ator T>i l in this paper, we obtain 

yh(TV)i - xh(rv)l = £• (3-82) 

Applying the decomposition formulae (|3-65p and ()3-72p . the gauge-transformation 
rules for the variable h(rvL) an d hrpvv) are given by 

yh(TVL) - xh(TVL) = £(L)h (3-83) 
yh(TW)i - xh(Tvv)i = i{v)i (3-84) 
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since we ignore the modes which belong to the kernel of A in this paper. Further, 
substituting Eq. (|3-82j) into Eq. (|3-80p . we obtain the gauge-transformation rule for 
h(TT)ij as follows 

(3-85) 



yh(TT)ij 



xh{TT)ij 



0. 



Thus, the gauge transformation rule for the variables h u , h(L) : 

h(TV)i-> (h{TVL) an d hi T vv)i)-> an d ^(TT)ij are summarized as : 



y">tt 



yh(VL) - 
yh(v)i 



X^{VL) 

- xh(y)i 



d t i {L) + 6 + A' 1 [2A (K ij Dj£ iL) ) + D k Ki (v)k 
d t i (y)l + 2K j l D j i {L) + 2K\i {v)3 



(3-86) 
,(3-87) 



-DiA 



-i 



2A (^^e ( L)) + ^ fc ^e ( y )fc 

yh(TV)l - xh(T V )l 
yh(TT)ij - xh(j>T)ij 

The equation (|3-90p is equivalent to the gauge-transformation rules 

yh(TVL) - xh(TVL) = £(Z,)i 

yh(Tvv)i - xh(rvv)i = £(v)J- 

Now, we construct gauge-invariant variables. First, Eq. (|3-9ip shows that the 
variable hrpT)ij * s itself gauge invariant. Therefore, we define the transverse-traceless 
gauge-invariant tensor as 



0. 



(3-88) 

(3-89) 
(3-90) 
(3-91) 

(3-92) 
(3-93) 



h, 



(3-94) 



\TT)iy 

Second, through Eqs. (|3-88p . (|3-92p . and (|3-93p . we can define the gauge-invariant 
variable for vector mode as 



Vi 



h 



(V)i 



d t h 



(TVV)i 



2K J 



(Djh( TVL ) + h( TVV )j) 



+D%A 



-i 



2 A [K ij D j h (TVL) )+D k Kh 



(TVV)k 



(3-95) 



Actually, it is straightforward to confirm that the variable i>i defined by Eq. (|3-95p 
is gauge invariant. Further, we can also confirm the divergenceless property of the 
variable z^, i.e., D l v, L = through the definition (|3 95|) and the formula (|3-70|) . 

Next, we consider the gauge invariant variables for scalar modes. Before doing 
this, we first construct the variable Xt in Eq. (|3-60p . Inspecting gauge-transformation 
rules (|3-87p . (|3-92p . and (|3-93p . we consider the combination 



X t := h 



(VL) 



d t h 



(TVL) 



A 



-i 



2 A {K l W j h {TVL) )+D k Kh 



(TVV)k 



(3-96) 



The variable X t defined by Eq. (|3-96p does satisfy the gauge-transformation rule 
(|3-6ip . We also define the variable Xi by 



Xi := h(r V)i = Dih^vL) + h(rvv)i- 



(3-97) 
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From this gauge-transformation rule (|3-90p . or equivalently Eqs. (|3-92p and (|3-93p . 
the gauge-transformation rule for the variable Xi defined by Eq. (|3-97p is given by 

y Xi - xX = (3-98) 

Now, we define the gauge invariant variables for the scalar mode. First, inspect- 
ing gauge-transformation rules (|3-61|) and (|3-86|) . we define the variable & by 

- 2<P := ha ~ 28 t X t . (3-99) 

Actually, we can easily confirm that the variable <P is gauge invariant. Second, in- 
specting gauge-transformation rules (|3-89p and (|3-98p . we define the gauge-invariant 
variable \P by 

- 2n& := h {L) - 2D i X i . (3-100) 

In terms of the gauge-invariant variables Xiji v ii ^> an d which are defined by 
Eq. (|3-94p . (|3-95]) . (|3-99p . and (|3-100p . respectively, and the gauge- variant variables 
X t and Xi which are defined by Eqs. (|3-96|) and (|3-97p . respectively, the original set 
{hu, hti, hij} of the components of the linear metric perturbation is given by 

h tt = -2$ + 2d t X t , (3-101) 

h ti = Vi + D % X t + d t X t + 2K\Xj , (3- 102) 

= -2& qij + xij + DiX,j + DjXi + 2KijX t . (3-103) 

On the other hand, we consider the decomposition formula (|2-18p . The components 
of the expression (|2-18p in the situation of this subsection with 

X a =:X t (dt) a + X t (dx l ) a (3-104) 

are given by 

ha = n tt + 28 t X t , (3-105) 

h ti = Hti + d t Xi + DiX t + 2K\Xj, (3-106) 

Kj = Hij + D t Xj + DjXi + 2KijX t . (3-107) 

Comparing Eqs. (|3-10ip - (|3-103|) and Eqs. (|3-105|) - (|3-107p . we easily see a natural 
choice of the components of the gauge-invariant part H a b and the components of the 
gauge-variant parts X a are given by 

Utt = -20, Uti = n, Hij = -2Vqij+ X ij, (3-108) 

and 

X t = X t , X % = Xi. (3-109) 

Of course, we may add the Killing vectors associated with the metric q a b to the 
definition of X a . The gauge-transformation rules (|3-6ip and (|3-98p support these 
identifications. Thus, Eqs. (|3-108|) and (|3-109p show that the linear-order metric 
perturbation h a b is also decomposed into gauge-invariant and gauge-variant parts as 
Eq. (|2-18p even in the case for the background spacetime with a = 1, (3 l = 0, but 
Kij y£ 0. These results are already reported in the previous letteil^ by the present 
author. 
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3.3. The case for arbitrary a, Pi, and K^j 

Now, we consider the most generic case of the metric (|3-ip where a / 1, /5 ! / 0, 
and Kij ^ 0. Considering the components of the metric perturbation h a b as Eq. (13- 2p . 
the gauge-transformation rules for these components are given by Eqs. (|3-3p - (|3-5|) . 
To do this, we first assume that the existence of the variables Xt and Xi whose 
gauge-transformation rules are given by 

y X t - x X t = (3-110) 

yXi - X X % = fc. (3-111) 

This assumption is confirmed through the construction of the gauge-invariant vari- 
ables for the linear-order metric perturbation below. Inspecting gauge-transformation 
rules (|3-3p - (|3-5p . we define the symmetric tensor field H a b whose components are 
given by 

Htt := ha + - (d t a + j3'Aa - PP 1 ^) X t 

+o?D i a - aP k D i p k - tf$Djo\ X h (3-112) 



3 , 

Hti ■= hi + - (Aa - pKy) X t 

+^ [^-o?K\ + P j p k K ki - ft D ia + aA/3 j ) Xj, (3-113) 

Hij := hj - ^K l3 X t + ^KijXk. (3-114) 

The gauge transformation rules (|3-3|H|3-5p and our assumptions (|3-110|) and (|3-111|) 
give the gauge-transformation rules of the components of H a b as follows: 

yHu ~ xH tt = 2d£ t , (3-115) 
y H ti - x H ti = dtti + (3-116) 
yHij - xHij = 2D( i £j). (3-117) 

Since the components Hu and H^ are regarded as a vector and a symmetric 
tensor on U(t), respectively, we may apply the decomposition reviewed in Appendix 
|B]to Hu and Hij: 

Hu = Dih^ VL ) + h(y)i, D % h(y^i = 0, (3-118) 
Hij = ^Qij h (L) + h (T)ij: Q l3 h(T)ij = 0, (3-119) 

h(T)ij = {Lh^TV))^ + h (TT)ij, D l h(T T )ij = 0, (3-120) 
h(TV)i = D ih(TVL) + h(TVV)h D % h( TV V)i = 0- (3-121) 

The gauge-transformation rules (|3-116p and (|3-117p gives the gauge-transformation 
rules for the variables h^ VL ^, hry^, fyn, ^(T)ij> h(TV)i ( or equivalently hrrvL) an d 
h(rw)i)i an d ^(TT)ij as i n the previous subsection. 
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First, we consider the gauge-transformation rule (|3T16p in terms of the decom- 
position (|3-118D : 

yHu - xH t i = A {yh(vL) ~ xh(VL)) + {yh(v)i ~ yh(v)i) = 3t& + A6-(3T22) 

Taking the divergence of this gauge-transformation rule and through the property 
D l h(y}i = 0, we obtain 

^ [yh ( vL) ~ xHvl)) = D%Si + (3-123) 

To evaluate the first term D l dt(,i in Eq. ()3-123[) . we note that dtq%j is given by 

d mj = -laKij + 2D (i (3 j} (3-124) 

from Eq. (|A-37p in Appendix [XJ Keep this equation in our mind, we consider the 
derivative Djdt^i as in the previous subsection. This is given by 

= d t DjCi - a {aK kj - A fc /y e fc 

-Dj {ptK kl - D {k ^) i k + A (oKy - A*/?;)) £*• (3-125) 
From Eq. (|3T25|) . we can easily derive D l dt£,i as 

D'dtti = d t D% - 2 (aK ij - A'ft')) D& 

-2A (aK H - A'ft)) 6 + D l (aK - Aft) (3-126) 
Further, we consider the decomposition of the component £j as Eq. (|3-72|) . i.e., 

6 = Ae ( L) + £(v)i> ^V)* = °- ( 3 ' 127 ) 

Through the decomposition ()3-127|) and the formula (|3-126p . the gauge-transformation 
rule (|3T23p is given by 



4 {yh {VL) - x h {V L) ~ it) = Ad£ [L) - 2 [aK^ - D^) D^ (v)l 

-2A (aK li - A' ft)) i {v)l 

+A (aK - Aft) i iv)l , (3-128) 

where we have used Eq. ()3-126|) twice. Ignoring the mode which belongs to the kernel 
of the derivative operator A, we obtain 

yh(VL) ~ X h (VL) =Ct + dtt,(L) 

+A- 1 [-2 (aK ij - A 4 ft)) D^i 
+ {A (aK- Aft) 

-2A (aK H - A'ft>) } C( V )l] ■ (3-129) 



24 



K. Nakamura 



Substituting Eq. (|3-129p into Eq. (|3-122p we obtain 

yh(V)i ~ yh(V)i = dtt( V )i ~ D.A- 1 [-2 (aKV - D^ {v)i 

+ 



{d 1 (aK - D { Pi 



-2 A ( aK u - D 



il (3 l) )}t { V)i] -(3-130) 



The gauge-transformation rules for /i^ and /i(T)«j are given from Eq. (|3-117p . 
Since we consider the decomposition (j3- 119|) . the gauge-transformation rule ()3-117|) 
is given by 

yHij - xHij = -qij (yh(L) - xh(L)) + {yh(T)ij - xh(T)ij) = 2D (i£j)- (3-131) 

Taking the trace of Eq. (|3-13ip , we obtain 

yh {L) - X h {L) =2D%. (3-132) 

The traceless part of Eq. (|3-13ip is given by 

yh(T)ij ~ x^{T)ij = ( L £)ij • (3-133) 

Note that the variable h/j>)ij is a l so decomposed as Eq. (|3-120p and the gauge- 
transformation rules for the variable fyr)^ is given by 

yh(T)ij - xh( T )ij = (L (yh( TV ) - xh^v)))^ + yh(TT)ij - xh(rr)ij 

= (Li\ r (3-134) 
Taking the divergence of Eq. (|3-134|) . we obtain 

V jl {yh (TV)l - x h {TV) i - 6) = 0. (3-135) 
Since we ignore the modes which belong to the kernel of T>i l , we obtain 

yh(TV)l - xh(T V )i = ii- (3-136) 
Through the decomposition formula ()3-12ip and (|3-127p . we easily derive 

yh(TVL) - xhpvL) = £(£)> (3-137) 
yh(TVV)i - xh(TVV)i = C(V)i) (3-138) 

where we ignore the mode which belong to the kernel of A. Substituting Eq. (|3-136p 
into (|3T34p . we obtain 

yh(TT)ij - xh(TT)ij = 0. (3-139) 

In summary, we have obtained the gauge-transformation rules for the variables 
Htt, h( VL ), /i(L), ^(T)y> fyrv)i; h( TVL ), h^yyy, and /i(TT)ij as follows: 

y H tt - x H tt = 2d t Z u (3-140) 
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yh(VL) - xh(yL) = 6 + 9t£(L) 



(3-141) 



yh(V)i - xh(y)i — d t £,(v)i 



DiA 



y h {L) - 
yh(T)ij 



- xh {T )ij 



2D%, 

= £(L), 

yh(TVV)i - xh(TVV)i = £(V)i 
yh(TT)ij - xh(TT)ij = 0. 



yh{Tv)i 

yh(TVL) 



xh(rv)i = 

- xh(TVL) 



2 (aKV - D^pti 
+ |d z (aif - D%) 

-2A (ai^-Z^)}^" 
-2 (aif^ - D^) D^ (y)fe 

-2D fc (aK lk - D«f3 k ^ } (3-142) 

(3-143) 
(3-144) 
(3-145) 
(3-146) 
(3-147) 
(3-148) 



Here, we note that the gauge transformation rule (|3-144|) coincides with the 
gauge transformation rule (|3-llip for the variable Xj. Then, we may identify the 



variable X; with h 



(TV)i- 



Xi 



h 



(TV)i- 



(3-149) 



Thus, we have confirmed the existence of the variable JQ. Next, we show the ex- 
istence of the variable Xt whose gauge-transformation rule is given by Eq. (|3-110p . 
To do this, we consider the gauge transformation rules f|3- 141 [) . (|3-146p . and (|3-147p . 
Inspecting these gauge transformation rules, we find the definition of Xt as 



h 



(VL) 



9th(TVL) 



2 (aK^ — D^ftn Djh 



J j n {TVV)i 



+ {d 1 (aK - D^i) - 2D { (aK li - D^ l (3 1 ^ } fr (7W )i] • (3-150) 

Actually, the gauge transformation rule for Xt defined by Eq. (|3-150p is given by 
Eq. (|3-11U|) . This is desired property for the variable Xt. Thus, we have confirm the 
existence of the variables Xt and Xi which was assumed in the definitions (|3-112p - 
(|3-114p of the components of the tensor field H^. 

Now, we construct gauge invariant variables for the linear-order metric pertur- 
bation. First, the gauge transformation rule (|3-148p shows that htj<T)ii ls gauge 
invariant by itself and we define the gauge-invariant transverse-traceless tensor by 



Xij '■— h(TT)ij- 



(3-151) 
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Inspecting the gauge-transformation rules (|3-142p and (|3-147p . we define the 
vector mode i>i by 

I/j := hfy)i - 9th(TVV)i 

+D..A- 1 [-2 (aK k i - £>(*0>>) i^tyzw)* 

+ [D l (aK - D k h) 

-2D k (aK lk - L>(^ fe )) } h (TVV)l ] . (3-152) 

Actually, we can easily confirm that the variable z/j is gauge invariant, i.e., yVi— x v % = 
0. Through the divergenceless property D t h^ v ^ i = for the variable hry^, we easily 
derive 

D*Ui = -D l d t h [TVV)t + (-2aK k * + 2D^ k (3^ D 3 h {TVV)k 

+ {-2D k [aK lk - D (l p k ^ + D l (aK - D k (3 k )j } h {TVV)l . (3-153) 

Further, through the formula (|3-126p for the variable hrrvv)i an d the divergenceless 
property of the variable h(j>w)ii we easily see the divergenceless property -DVj = 0. 

Next, we consider the scalar modes. First, inspecting gauge-transformation rules 
(|3-110|) and (|3-115p . we define the scalar variable <P by 

- 2<2> := H u - 2d t X t . (3-154) 

Actually, we can easily confirm that this variable <£> is gauge invariant. Inspecting 
the gauge-transformation rules ()3-llip and (|3-132p , we define another gauge-invariant 
variable \P by 

- 2n& := h {L) - 2D i X i . (3-155) 

We can easily confirm the gauge invariance of the variable & through the gauge- 
transformation rules (|3-111|) and (|3-132p . 

In summary, we have defined gauge invariant variables as follows: 



-2<2> 

V; 



H u - 2dtX t , (3-156) 
h {L) - 2D l X u (3-157) 
h(V)i - dth(TVV)i 



+DiA- x [-2 (aK k 3 - D&pti) D 3 h {TVV)k 
+ {D l (aK-D k [3 k 
-2D h ( aK lk - 



D {l f3 k) )}h (TVV)l ] , (3-158) 



Xij ■= h (TT)ij- (3-159) 

In terms of these gauge- invariant variables and the variables Xt and Xi, which are 
defined by Eqs. (|3-15Up and (|3-149|) . respectively, the original components {h tt , hu, 
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h^} of the metric perturbation h a b is given by 

ha = -2$ + 2d t X t - - (d t a + pD % a - p?K iS ) X t 
-1 (j?pl>pK kj _ pi dta + aqiidtpj 

+a 2 £>*a - aP h D i p h - frfiWjoJ X h (3-160) 

h ti =Vi + D t X t + d t Xi - - (D.a - pKy) X t 

a 

~ (~a 2 K\ + P(3 k K ki - p3 Dia + a D^ Xj, (3-161) 

^ = -2W qij + Xi j + DiXj + DjXi + ^KijX t - ^(3 k KijX k . (3-162) 

On the other hand, the component representations of the decomposition formula 
(12481) with 



X a =: X t (dt) a + X,{dx l ) a (3-163) 



are given by 



U tt + 2d t X t - - (d t a + FDta - /3 fc /3^j) X t 



+a 2 Aa - a(3 k D l (3 k - FfiDjCt ) X t , (3-164) 



hi = Hti + d t X % + DiX t - - (A« - P' J Kij) X t 

-- [~a 2 K\ + P j (3 k K ki - p j Di,a + aA/3 J ) Xj, (3-165) 

hij = Hij + DiXj + DjXi + ^KijX t - ?-p k KijX k . (3-166) 

Comparing Eqs. (j3- 164[ >— (j3- 166p with Eqs. (|3-160p - (|3-162p . we may identify the com- 
ponents of the gauge- invariant variables H a b so that 

U U :=-2<P, n ti :=Ui, Hij := -2& qij + X ij (3-167) 

and the components of the gauge-variant variables X a so that 

X t := X t , X % := Xi. (3-168) 

Thus, the decomposition formula f)2-18[) is correct for the linear-order perturbation 
on a generic background spacetime, if we assume the existence of two Green function 
of the derivative operators A := D l Di and T> 10 which is defined by Eq. (|3-26|) . In 
other words, in the above proof, we ignore the modes which belong to the kernel of 
these derivative operators A and To take these modes into account, the different 
treatments are necessary. 
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§4. Comparison with the FRW background case 

In this section, we consider the comparison with the case where the back- 
ground spacetime Mq is a homogeneous and isotropic universe which is discussed in 
KN2007J-^ This case corresponds to the case a = 1, j3 l = and Kij = —Hqij, where 
H = dta/a and a is the scale factor of the universe. 

In the paper KN2007, 7 ' we consider the decomposition of the components hu 
and of the metric perturbation h a b as 



hi 


= ^>ih(vh) + h(y)i, D l h(y)i = 0, 


(4-1) 


hij 


= a 2 ^(L)7ij + a 2 h(T)ij, 7 u ^(T)ij = 0, 


(4-2) 


h(T)ij 


= (^DiDj - -^HjAj h(TL) + 2 D(ih(TV)j) + h(TT)ij, 


(4-3) 


D l h(r V )i 


= 0, D\ TT)i j = 0, 


(4-4) 



where = a 2 ^fij, 7™ is the metric on a maximally symmetric space, Di is the covari- 
ant derivative associated with the metric jij, and A := D l Di. This decomposition 
is slightly different from the decomposition (|3-62p - (j3-65p with the definition ()3-60p 
of the variable H a b- Furthermore, as noted in KN2007P there should exist Green 
functions of the derivative operators A, A + 2K, and A + 3K to guarantee the one 
to one correspondence of the set {h tt , h u , h^} and {{h tt , h (TL)}, {fyv)i> 

tyxv)t}j h(TT)ij}> where K is the curvature constant on the maximally symmetric 
space. The special modes which belong to the kernel of the derivative operators A, 
A + (n — 1)K, and A + nK were not included in the consideration of the paper 
KN2OO7. On the other hand, in this paper, we ignore the modes which belong to 
the kernel of the derivative operator A and T) 13 . In this section, we briefly discuss 
these correspondence. 

First, we note that the decomposition (|4-ip of the component ha of the met- 
ric perturbation h^b is equivalent to (|3-62p . Although the tiny difference between 
Eq. ()4-ip and (|3-62p is in the definition of the covariant derivatives Di (associated 
with the metric q^) and Di (associated with the metric ^ := (1/a 2 )%.,■), we may 
say that fyyz,) and h(yy in Eq. (|4-ip are identical with h(y L ^ and h(yy in Eq. (|3-62p . 
respectively. 

We also note that the trace parts of these two decompositions are almost equiva- 
lent. Actually, since the extrinsic curvature Kij on the background U is proportional 
to the intrinsic metric q^ in this case, we easily see that 

'/"/,,, = q lj Hij + 2q i iK ij X t 

= h (L) - 2nH (h {VL) - d t h {TVL) + 2Hh {TVL) ) , (4-5) 

where we used Eq. (|3-96p . On the other hand, the trance part of hij given by Eq. (|4-2p 
is nh^y Thus, the variable h^ in Eq. (|4-2p corresponds to the variables in this 
paper as 

h{L) = -h( L ) ~ 2i? [h(VL) ~ dth(rvL) + 2Hh {TVL) ) . (4-6) 
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Since the extrinsic curvature in Eq. ()3-60p is proportional to the intrinsic 
metric q^ in this case, the main difference between decompositions (|3-63j) (|3-65|) 
and (|4-2|) - (|4-4p are in the traceless part. The traceless part hrp)ij in Eq. (|4-3p is also 
given by 

h(T)ij = f A-Dj - n lij ^) H{TL) + 2 ^i^ TV )i) + ~ h (TT)ij 

-^q ij a 2 D k (^D k h(r L ) + h {TV) ^j + a 2 h {TT)ij , (4-7) 



where we used 

= D k h {TV)k = a 2 D k h {TV)k . (4-8) 
Comparing (|3-64p . we obtain the correspondence of the variables 

h(TV)i = 2 a2 A^(TL) + a2 h(TV)ji ^(TT)ij = a2 h(TT)ij, (4-9) 

h(TVL) = ^^(tl), fyTVV)i = a2 h(rv)i- (4-10) 

Therefore, in the case of the homogeneous isotropic universe, the decomposition 
(JM^ - fpHlSj) is equivalent to the decomposition (pFTjlHTI^j) . 

However, in the case of the generic background spacetime, the decomposition 
(|4-1|) - (|4-4|) is ill-defined. Actually, if we regard that the decomposition (|4-ip - (|4-4|) 
is that for the generic background spacetime, we cannot separate /i(TL) an d h(TV)j 
due to the non-trivial curvature terms of the background Mq as pointed out by 
DeserPS These curvature terms come from the commutation relation between the 
covariant derivative -Dj and the derivative operator D lJ . This is why we apply the 
decomposition (p^U2l - (l3^1) instead of (|4"TD -(I¥4D . 

Finally, we consider the correspondence of the special modes which we ignore 
in this paper and KN2007.^ Trivially, the above operator A := D % Di corresponds 
to the Laplacian A in this paper. The above derivative operator A + (n — 1)K 
corresponds to the derivative operator T> lJ . In the case of the maximally symmetric 
n-space, the Riemann curvature and Ricci curvature are given by 

(n) Ri jk i = 2Kq k[iqj]l = 2Kq k[i q j]b R lk = q jl{n) R ijk i = (n - l)Kq ik . (4-11) 

In this case, the derivative operator T> %3 defined by Eq. f|3-26j) is given by 

V ij = q ij (A + (n- l)K) + M - ^ J D { D j . (4-12) 

When the operator T> %1 acts on an arbitrary transverse vector field Vi {D l Vi = 0), we 
easily see that 

V ij Vj = (A + (n- l)K) v\ (4-13) 
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Finally, we point out that the above derivative operator A + nK appears in the case 
where the derivative operator T>- 1 acts on the gradient Dif of an arbitrary scalar 
function /. Actually, we easily see that 



Tfl'Dif = 2- 



n 



D j A + — - — B? l Di 
n — 1 



/• (4-14) 



In the case of maximally symmetric n-space, curvature tensors are given by Eqs. (|4-lip 
and the derivative operator T> 3l Di is given by 

V jl Dif = 2^—^-D j (A + nK)f. (4-15) 

n 

When we solve the equation 

Tfi l Dif = gi, (4-16) 

we have to use the Green function A and A + nK. These are the reason for the fact 
that the Green functions A~ 1 , (A + (n — 1)K)~ , and (A + nK)~ l were necessary 
to guarantee the one-to-one correspondence between the components {hu, hij} and 

{h(vL), h(y)i, h(L), h(TL) : h(TV)i> h(TT)ij} in Eqs. (gd])-(|131). In other words, we 
may say that the special modes belong to the kernel of the derivative operators A 
and T> %3 which are ignored in this paper are equivalent to the special modes which 
belong to the kernel of the derivative operators A, A + {n — 1)K, and A + nK which 
are ignored in the paper KN2007.® 



§5. Summary and discussions 



In summary, after reviewing the general framework of the higher-order gauge- 
invariant perturbation theory in general relativity, we prove Conjecture 12.11 for 
generic background spacetime which admits ADM decomposition. In this proof, 
we assumed the existence of Green functions of the elliptic derivative operators A 
and T> %0 . Roughly speaking, Conjecture 12.11 states that we know the procedure to 
decompose the linear-order metric perturbation h a b into its gauge-invariant part Tl a b 
and gauge- variant part X a . In the cosmological perturbation case, this conjecture 
is confirmed and the second-order cosmological perturbation theory was developed 
in our series of papers.^'^'^ 1 However, as reviewed in SJ21 Conjecture 12.11 was the 
only non-trivial part when we consider the general framework of gauge-invariant 
perturbation theory on generic background spacetimes. Although there may exist 
many approaches to prove Conjecture 12.11 in this paper, we just proposed a proof 
for generic background spacetimes. 

As noted above, in our proof, we assume the existence of the Green functions for 
the elliptic derivative operators A and T) 13 . This assumption implies that we have 
ignored the modes which belong to the kernel of these derivative operators. Within 
the arguments in this paper, there is no information for the treatment of these mode. 
To discuss these modes, different treatments of perturbations are necessary. We call 
this problem as zero-mode 'problem. The situation is similar to the cosmological 
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perturbation case as noted in £§] and zero-mode problem exists even in the cosmo- 
logical perturbation case. In the cosmological perturbation case, zero-mode means 
the modes which belong to the kernel of the derivative operator A, A + (n — 1)K, 
and A + nK, where K is the curvature constant of the maximally symmetric space 
in cosmology and n is the dimension of this maximally symmetric space. 

This zero-mode problem in cosmological perturbations also corresponds to the 
I = and I = 1 mode problem in perturbation theory on spherically symmetric 
background spacetimes. In the perturbation theory on spherically symmetric back- 
ground spacetimes, we consider the similar decomposition to Eqs. (|4-2p - (|4-4p and 
the indices in these equations correspond to the indices of the components of 

a tensor field on S 2 . Since S 2 is a 2-dimensional maximally symmetric space with 
the positive curvature, we may regard n = 2 and K = 1. Then, the above three 
derivative operators are given by A, A + 1, and A + 2. Since the eigenvalue of the 
Laplacian A on S 2 is given by A = —1(1 + 1), we may say that the modes with I = 
and I = 1 belong to the kernel of the derivative operator A, A + (n — 1)K, and 
A + nK. Therefore, we may say that the problem concerning about the modes with 
I = and I = 1 in the perturbations on spherically symmetric background spacetime 
is the same problem as the zero-mode problem mentioned above. 

Thus, the arguments in this paper shows that zero-mode problem generally ap- 
pears in many perturbation theories in general relativity we have seen that the 
appearance of this zero-mode problem from general point of view. To resolve this 
zero-mode problem, carefully discussions on domains of functions for perturbations 
will be necessary. We leave this zero-mode problem as a future work. 

Although we should take care of the zero-mode problem, we have almost com- 
pleted the general framework of the higher-order gauge-invariant perturbation theory 
in general relativity. The proof of Conjecture 12. ll shown in this paper gives rise to the 
possibility of the application of our general framework for the higher-order gauge- 
invariant perturbation theory not only to cosmological perturbations^'^'^ 1 but 
also to perturbations of black hole spacetimes or perturbations of general relativis- 
tic stars. Therefore, we may say that the wide applications of our gauge-invariant 
perturbation theory are opened due to the discussions in this paper. We also leave 
these development of gauge-invariant perturbation theories for these background 
spacetimes as future works. 
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Appendix A 

ADM decomposition 

Here, we briefly review the ADM decomposition.^'^ 
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We consider the n + 1-dimensional spacetime (Ai, g a b)- The topology of Ai is 
given by Ai = M 1 x E, where E is the n-dimensional manifold. This means that the 
entire Ai is foliated by the one-parameter family of the manifolds Et where t is the 
parameter along M 1 in M 1 x E. Here, we note that it is not necessary to impose that 
the entire spacetime Ai is decomposed into M 1 x E in the global sense. However, 
in this section, we impose that there exists a one-parameter family of n-dimensional 
submanifolds Et in a n + 1-dimensional manifold Ai , for simplicity. 

In general relativity, t is regarded as the time function on Ai and the decompo- 
sition of Ai into M 1 x E is regarded as the n + 1-decomposition of the spacetime Ai 
into the space E and time M 1 . Prom the view point of this decomposition, we can 
describe the Ai by the time-evolution of E, i.e., the geometry of Ai is described in 
terms of the geometry of E which is embedded in Ai. Any geometrical quantities 
on Ai is given in terms of the geometry of E and its "time-evolution" . 

First, we consider the n + 1-decomposition of the metric g ab on Ai. Let t a be a 
vector field on Ai satisfying 

t a V a t = 1, (AI) 

i.e., 




(A-2) 



Let us call the direction along which the function t increases as the future direc- 
tion and t a defined by (|A-1[) is called future-directed, the direction along which the 
function t decrease as the past direction. 

On the other hand, let us denote the unit normal to Et by n a , which is hyper- 
surface orthogonal. The normalization condition for n a is given by 

n a n a = -1. (A-3) 

The metric g a b on Ai induces a metric q a b on E t . On each E t , q a b is given by 

Qab ■= 9ab + n a n b . (A-4) 

The overall signature of n a is chosen so that n a is future-directed. Then there is a 
positive function a so that 

n a V a t=:-. (A-5) 

a 

This positive function a is called the lapse function with respect to t a . Further, we 
impose the n a is the hypersurface orthogonal to the hypersurfaces Ef, which implies 
that n a oc V a t. Due to the normalization condition (|A-3p . we easily see that 

n a = -aV a t = -a (dt) a . (A-6) 

From Eq. ()A-6p . we decompose the vector field t a into its normal and tangential 
parts to Et 



a = -t a n a , P a := q ab t b , 
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where (3 a is called the shift vector with respect to t a . Equivalently, the vector field 
t a is decomposed as 

t a = t b 5i = p a + an a . (A-8) 
From Eq. (|A-1|) and the definition (|A-5p of the lapse function a, we easily see that 

1 = t a V a t = an a V a t + /3°V t = 1 + /3°V t, (A-9) 

which yield 

/3 a V a t = 0, (A-10) 

i.e., f3 a has no component along (d/dt) a . Further, due to the normalization (|A-3j) of 
the vector n a , the lapse function a is also given by 

a 2 = \ • (A-ll) 

The decomposition (|A-8[) of the vector field i a = (d/dt) a is also yields 

n° = i [t° - )9 a l . (A- 12) 

a 

Introducing the spatial coordinate so that (t,x l ) is the spacetime coordinate, 
i.e., the coordinate basis of the tangent space on the spacetime is the set 

{(*)„(*").}, {(!)'. (£)'}. < A13 > 

/3 a is given by 



because of Eq. (|A-10j) . Then, (|A-12p is given in terms of the coordinate system 

1 



n = — 

a 



at/ p \dx i 



(A-15) 



From 

V a t = 5 a b V b t = (g a c - n a n c ) V c tq a b V b t -n a -. (A- 16) 

a 

and (|A-6[) . we can also see that 

g a b V b t = 0. (A-17) 

This implies 

q ab V b t = q ba V b t = 0, (A-18) 

which yields q ab as no component along (d/dt) a . This means that the induced inverse 
metric q ab has the following component representation 
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Together with the component representation (|A-12p and (|A-19p . we obtain the space- 
time inverse metric g ab in terms of the coordinate system (t,x i : 



g ab = _ n a n b + g a6 j ( A . 2() ) 



a 2 \\dt) 13 \8x'J j \ \8t) P \dxi) 
Through the relation 

9ab9 bc = S a c , (A-22) 

the straight forward calculation leads the coordinate representation of the metric 
gab, which is given by 

g ab = -a 2 {dt) a (dt) b + qijidx 1 + Pd£) a {dxi + Fdt) b , (A-23) 

where (dt, dx l ) is the coordinate basis on Ai (more precisely, on an open set U C Ai), 
and qij is the inverse matrix of q %3 in Eq. (|A-2ip . i.e., 

qij gjk = . q .k = S .k (A . 24) 

and (5j k is the n-dimensional Kronecker's delta. In terms of the coordinate basis 
(|A-23p . the unit normal vector n a is given by 

n a = -a(dt) a , (A-25) 

Here, we summarize the components of the spacetime metric g ab and inverse 
metric g ab on Ai as follows: 

g tt = -a 2 + qijP l /3 J , gu = 9u = Qijft = A, 9%j = qij, (A-27) 

9* = -\, g u = g lt = \p\ .'/'•' '/'•' -, '>' '>■'■ (A-28) 

or or ot z 

Next, we consider the connection between the covariant derivative V a on (Ai,g ab ) 
and the covariant derivative D a on (E,q ab ). This correspondence is given from the 
Christoffel symbol associated with the metric g ab in the coordinate system (|A-23p : 

1% = -d t a + -pD ia + ^ k p l {d t qij - 2D {1 (3 J} } , (A-29) 
T \ t = I Di a + ^2 P j {dtqij - 2D {i /3;) I . ( A - 30 ) 
4 = ^2 {^-2^)}, (A-31) 



a 
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+aD l a - l3 k D l (3 k - -ppDja, (A-32) 



1 



a 

Jk fa „ on fl 1 ^ oi pk 



r h = 2^ - 2%^)} - ^2 ?P i 9 ^ ~ 2D (fc /5tf)} 

-pDja + Djp, (A-33) 



a 



= {%fej - 2%^)} + (A-34) 
where ^-T^ is the Christoffel symbol associated with the metric g^: 

&>rj k = (d jqik + d m - d iqjk ) . (A-35) 

It is also convenient to introduce the extrinsic curvature 

K ab = -q a c q b d V c n d . (A-36) 

In the coordinate system, on which the metric is given by (|A-23|) . this extrinsic 
curvature is given by 

d 



Ka = — 
3 2a 



(A-37) 



Through this component of the extrinsic curvature Kij, the above components 
(|A-29|H|A-34|) of the Christoffel symbols T 6 a c are given by 

r t \ = -d t a + VAo - -P k ^Kij, (A-38) 
a a a 

Tl = -Did - K i:i , (A-39) 
a a 

4 = —Ka, (A-40) 



a a 

+aD' l a - P k D i p k - -ppDjOi, (A-41) 
a 

rL = -aK\ + ~pp k K ki - -PDiCt + Djfi, (A-42) 
J J a J a 

r; k = -^K k3 + ^r; k . (a-43) 

J a 

Appendix B 

Covariant orthogonal decomposition of symmetric tensors 

Since the each order metric perturbation is regarded as a symmetric tensor on the 
background spacetime (Mo, gab) through an appropriate gauge choice, the covariant 
decomposition of symmetric tensors is useful and actually used in the main text. 
Here, we review of the covariant decomposition of symmetric tensors of the second 
rank on an curved Riemannian manifold based on the work by YorkPSJ 
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In the generic curved Riemannian space (U, q a b) (dim 17 = n), one can decom- 
pose an arbitrary vector or one-form into its transverse and longitudinal parts as 

A a = Ai(dx% = (DiA [L) + A {v)i ) (dx%, D i A(y) i = 0, (B-l) 

where Di is the covariant derivative associated with the metric q a b = qij{dx t ) a {dx :! )b- 
-A(L) is called the longitudinal part or the scalar part and Ary)i is called the transverse 
part or vector part of the vector field A a on (£, q a b), respectively. 

Moreover, this decomposition is not only covariant with respect to arbitrary 
coordinate transformations, it is also orthogonal in the natural global scalar product. 
To clarify this orthogonality, York® introduced the inner product for the vector 
fields on U. This is, for any two vectors V a and W a , we have 

e q V a W b q ab , (B-2) 

where e q denotes the volume element which makes the integral invariant and the 
integration extends over the entire manifold (£, q a b)- In terms of this inner product, 
the orthogonality of the vector fields V a = D a Vii\ := q ab D b V^ and W a = q ab V^ V )b 
with D a V(v) a = is given by 

/ e q D a V {L) V (v)b q ab = [ s a V (L) V (v)b q ab - [ e q V [L) D a V (v)b q ab , (B-3) 

where s a is the volume element of the (n — l)-dimensional boundary dU of S. 
Since the second term of Eq. (jB-3[) vanishes due to the condition D a Vry\ a = 0, the 
inner product (V, W) vanishes if V(n and V(v)b satisfy some appropriate boundary 
conditions at the boundary dU of U so that the first term of Eq. (|B-ip vanishes. In 
this sense, the scalar part (the first term in Eq. (|B-ip ) and the vector part (the second 
term in Eq. (|BTj) ) orthogonal to each other. Geometrically, the decomposition of 1- 
forms, and more generally p- forms, leads via de Rham's theorem to a characterization 
of topological invariants of U (i.e., Betti Numbers).® 

In this appendix, it is assumed that the n-dimensional space U is closed (com- 
pact manifolds without boundary) following York's discussions. The choice of closed 
spaces is made for mathematical convenience but the decomposition discussed here 
is also valid for any other n-dimensional spaces U with the boundary dS with some 
appropriate boundary conditions at dE. Through this assumption, in this appendix, 
we consider the TT-decomposition (transverse traceless decomposition) of a symmet- 
ric tensor t/) ab on S, which is defined by 

r b = ^ b T + ^t+^T b n (B-4) 

where the longitudinal part is 

$f ■= D a W b + D b W a - -q ab D c W c =: (LW) ab (B-5) 



and the trace part is 



ijjf r :=-^q ab , ^:=q cd i> cd . (B-6) 
n 
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Let us suppose that both an arbitrary symmetric tensor field ifj ab and the metric 
q a b are C°° tensor fields on E. First, we define tp^ in accordance with Eq. (|B-4p by 

■= - -^ g ab - (LW) ab . (B-7) 
n 

We note that the tensor ip^ T is traceless, i.e., 

Qab^TT = (B-8) 

by its construction (|B-7p , Further, we require the transversality on the tensor field 

D b ^f T = 0. (B-9) 



Equation (|B-9p leads to a covariant equation of the vector field W a in Eq. (|B-7p as 

D a {LW) ab = D a (^ ab - ^q ab ^j • (B-10) 
The explicit expression of ()B-10p is given by 

V bc W c = D a {^ ab - ^q ab ^j , (B-11) 

where the derivative operator T> bc is defined by 

V bc := q bc A + ^1 - ^) D b D c + R bc , A := D a D a . (B-12) 

The basic properties of Eq. (|B-lip are also discussed by YorkP^ The operator 
T> ab defined by Eq. (|B-12p is linear and second order by its definition. As discussed 
by York, this operator is strongly elliptic, negative-definite, self-adjoint, and its "har- 
monic" functions are always orthogonal to the source (right-hand side) in Eq. (|B-lip . 
Here, "harmonic" functions of T> ab means functions which belong the kernel of the 
operator T> ab . Moreover, he showed that Eq. (|B-lip will always possess solutions 
W a which is unique up to conformal Killing vectors. Due to these situation, in this 
paper, we assume that the Green function (P -1 ) a & defined by 

(V-'UV^ = V bc (V-\ b = C (B-13) 

exists through appropriate boundary conditions at the boundary dS of U. Although 
York's discussions are for the case of the closed space U, we review his discussions 
here. In this review, we explicitly write the boundary terms which are neglected 
by the closed boundary condition to keep the extendibility to non-closed £ case of 
discussions in our mind. 

The ellipticity of an operator depends only upon its principal part, i.e., the 
highest derivatives acting on the unknown quantities which it contains. To see 
the ellipticity of an operator, we consider the replacement of the each derivative 
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operator D a occurring in its principal part by an arbitrary vector V a . Through this 
replacement, the principal part of the operator defines a linear transformation a v . 
The operator is said to be elliptic if o~ v is an isomorphism. 2 ?.' In the present case, 



ab 



ybya + q aby c yc^ 



(B44) 



Here, a v operates on vector X a and defines a vector-space isomorphism when the 
determinant of a v is non- vanishing for all non- vanishing V a . The fact that det a v 7^ 
here is verified, for example, by choosing V a = (d/dx ll ) a in a local Cartesian frame 
The operator is said to be strongly elliptic if all the eigenvalues of a v are 
nonvanishing and have the same sign. This is easily checked and T> ab is strongly 
elliptic. 

To show that T> ab is negative definite, we consider the inner product (|B-2j) of the 
vector field VW a := V ab W b and W a : 

f^ q q ab W a V bc W c = f^e q q ab W a (D c {LW) bc ) 

= J^e q Id c (w b (LW) bc ) - X -{LW) bc {LW) 
s c W b {LW) bc - \ ! e q (LW) bc (LW) 



be 



or 



be 



(B-15) 



where we use the fact that the tensor (LW) bc is symmetric and traceless. Eq. (|B-b5|) 
shows that the operator T> ab has the negative eigenvalues in the case where the 
first term (boundary term) in Eq. (|B-15|) is neg lected, unless {LW) bc = 0. The 
self-adjointness of the operator T> ab is follows from a similar argument in which one 
integrates by parts twice: 



e q q ab V a (VW) b = f e q q ab V a (D c {LW) bc 



D c 


[q ab V a (LW) bc ) 


D c 


(q ab V a (LW) bc ) 


D c 


(q ab V a (LW) bc ) 


D c 


(q ah V a {LW) hc ) 


D c 


[q ab V a {LW) bc ) 



be 



V c V b {LW) 



\ (LV) bc (LW) bc 



(LV) bc D c W b 



D c ((LV) bc W b ) +W b V c (LV) 



D c ((LV) bc W b ) +W b V bc V c 



V b (LW) bc - (LV) oc W b 



be * 



+ 



E 



e q W b V bc V c 



(B-16) 



for any vectors V and W, where we use the fact that the tensor (LW) ab and (LV) ab 
are symmetric and traceless. Eq. (|B-16p shows that the operator T> ab is self-adjoint 
if the first term (boundary term) in Eq. (|B-16p is neglected. 
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When we can neglect the boundary terms in Eq. (|BT5p . the right-hand side 
of (|Bd5|) can vanish only if (LW) ab = 0. This means either W a = or W a is a 
conformal Killing vector (or Killing vector) of the metric q a i>. The condition for a 
conformal Killing vector is, of course, not satisfied for an arbitrary metric but this 
is given by 

£wQab = Mob ( B-17 ) 

for some scalar function A, where £\y denotes the Lie derivative along W. Taking 
the trace of both sides, we find 

A = -V C W C . (B-18) 

n 

Therefore, W a is a conformal Killing vector if and only if 

V a W b + V b W a - -q ab V c W c = (LW) ab = 0. (B-19) 
3 

It follows that the only nontrivial solutions of T> ab W b = are conformal Killing 
vectors if they exist. Hence the nontrivial "harmonic" functions of T> ab are conformal 
Killing vectors. We shall now show that even if these "harmonic" solutions exist, 
they are always orthogonal to the right-hand side of (|B-10p and, hence, can cause no 
difficulties in solving equation (|B-10p by an eigen function expansion. 

Denote the conformal Killing vectors by W a = C a , where by definition (LC) ab = 
0. Form the scalar product of the right-hand side of (|BT0p with C and integrate by 
parts to find 



E 

e, 

E 



A, ( r A „. ( r b - C c ] - (V - ~<?°V 1 DuC, 



= s b Lf* - V 6 V>) C a - l -J^ e q - (LC) ab = 0, (B-20) 

where we use the fact that ip ab — ^q ab ifj is symmetric and traceless and we also neglect 
the boundary term. Hence the source in Eq. (|B- llj) is in the domain of (T)~ 1 ) ab 
and (T>~ 1 ) ab gives the solution to Eq. (jBTip even in the presence of conformal 
symmetries. 

These results also show that the solution to Eq. (jBTip must be unique up to 
conformal Killing vector fields. Since only [LW) ab enters in the definition (|B-7p of 
t^TT' conformal Killing vectors cannot affect ip^ T . 

The orthogonality of ip^ T , (LW) ab , and ^ipq ab is easily demonstrated. We see 
readily that ^ipq ab is pointwise orthogonal to (LW) ab and to ip^ T , as (LW) ab and 
iPtt are both trace-free. To show that ip^T an< ^ {LV) ab are orthogonal for any vector 
V and any TT tensor, we have only to show that 



'E 



e q qacqbd(LW) ab ip c T d T = J e q 2D a W b ^f T 
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= J e q (p a (2W h ^f T ) - 2W b D a il#r) 

= J s a {2W b ^f T ) - J e q (2W b D a ^f T ) = 0,(B-21) 

where we use the fact that the tensor i/j^t ^ s symmetric, traceless, and transverse 
(|B-9p . We also neglect the boundary term in Eq. (|B-2ip . Thus, we conclude that 
the decomposition defined by (|B-7p exists, is unique, and is orthogonal. 

One can further decompose the vector W a uniquely into its transverse and lon- 
gitudinal parts with respect to the metric q a b- This splitting is orthogonal, as in 

Eq. dnuD- 

Since the above discussions are for closed spaces E, careful discussions on the 
boundary terms which are neglected in the closed £ is necessary if we extend the 
above arguments to non-closed S case. However, we do not go into these detailed 
issues. Instead, in the main text, we assume that the existence of the Green function 
of the derivative operator T> ab and use the transverse-traceless decomposition for an 
arbitrary symmetric tensor on E discussed here. 
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